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Abstract. NGC 1866 is a young, rich star cluster in the Large Magellanic Cloud. Since the cluster is very well 
populated both in the main sequence and post main sequence stages, thus providing us with a statistically com- 
plete sample of objects throughout the various evolutionary phases of intermediate mass stars, it represents a good 
laboratory for testing stellar evolutionary models. More precisely, NGC 1866 can be used to discriminate among 
classical stellar models, in which the extension of the convective regions is fixed by the classical Schwarzschild 
criterion, from models with overshooting, in which an "extra-mixing" is consi dered to take place beyond the 
classical limit of the convective zone. Addressing this subject in a recent work, Testa et al. (1999) reached the 
conclusion that the classical scheme for the treatment of convection represents a good and sufficient approximation 
for convective interiors . Using their own data, we repeat here the analysis. First we revise the procedure followed 
by Testa et al. (1999) to correct the data for completeness, second we calculate new stellar models with up- 
dated physical input for both evolutionary schemes, finally we present many simulations of the colour-magnitude 
diagrams and luminosity functions of the cluster using the ratio of the integrated luminosity function of main 
sequence stars to the number of giants as the normalization factor of the simulations. We also take into account 
several possible physical agents that could alter the color-magnitude diagram and the luminosity function: they 
are unresolved binary stars, dispersion in the age, stochastic effects in the initial mass function. Their effect is 
analyzed separately, with the conclusion that binary stars have the largest impact. The main result of this study 
is that the convective overshoot hypothesis (together with a suitable percentage of unresolved binaries) is really 
needed to fully match the whole pattern of data. The main drawback of the classical models is that they cannot 
reproduce the correct ratio of main sequence to post-main sequence stars. 
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1. Introduction 

In the context of stellar theories, it has long been acknowl- 
edged that the penetration of convective elements beyond 
the classical limit set by the Schwarzschild criterion could 
produce non-negligible effects on the stellar structure. 

The young star cluster NGC 1866 in the Large 
Magellanic Cloud (LMC) is considered as an ideal labora- 
tory for testing stellar models, especially concerning the 
extension of convection in the interiors of real stars, and 
discriminating among classical and so-called overshoot- 
ing models. This cluster, well populated both in main se- 
quence and giants stars, represents a statistically complete 
sample of objects, contrary to the young clusters of the 
Milky Way (that are relatively poor, especially in giant 
stars). Furthermore, NGC 1866 is young enough (~100 
Myr) to possess stars which develop a convective core in 
the phase of central H-burning (Mro =4 — 5 M Q ), thus 
providing us the way to test the extension of overshooting. 
The theory of stellar structure predicts that convective 



cores in main sequence stars set in starting from initial 
masses greater than M ~ 1.1 M Q . For these reasons NGC 
1866 has been the subject of several analyses, aiming at 
testing how far convective elements overshoot from the 
core into the surrounding stable layers. 



Becker & Mathews (1983), comparing observational 
and synthetic colour-magnitude diagrams (CMD) of NGC 
1866, obtained the best fit to the data simulating a clus- 
ter of about 86 Myr: the general features of the CMD 
were well-fitted on the whole, but the authors argued that 
the classical stellar models they used predicted a number 
of red giants larger than observed, and a smaller num- 
ber of main sequence stars in turn as compared to the 
observations. The same authors suggested that a more 
careful treatment of core convection (i.e. larger convective 
cores) could remove the discrepancy. A number of subse- 
quent studies confirmed this idea and demonstrated that 
the inclusion of overshooting in the description of convec- 
tive motions could r eproduce the correct ratio ./Vpms/^Vms- 
|Chiosi et al. (1989a ), in particular, clearly showed that the 
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overshooting scheme, by reducing the ratio tn e /tn, offers a 
good and simple solution to the problem. This conclusion 
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was also reinforced by [Lattanzio ct al. (1991 ), who ana- 
lyzed the CMD of NGC 1866 by means of overshooting 
and classical models Q 

Over the years, sophisticated formulations of convec- 
tive overshooting have been elaborated that are based on 



to recall here the so-called mass discrepancy of Cepheid 



turbulence theories (Cloutman & Whitaker, 1980; Xiong 



1980) md fluid hydro-dynamics ( Canuto fc Mazzitclli 
199lF1( ianuto et al., 199"6| ; |Unno fc Kondo, 1989| ). However' 



the ballistic approach at the problem proposed by pressan 
et al 



1981 



turns out to be fully adequate and it has 
been proved to best reproduce the nu merical resu lts of 
laboratory fluid-dyn amics simulations ( Zahn, 1991 ). The 
Brcssan et al. ( 1981 ) algorithm adopts a non-local treat- 



ment of convection in the context of the mixing- length the- 
ory (MLT) by Bohm-Vitcnsc (1958 ) : it looks for the layers 
where the velocity of convective elements (accelerated by 
the buoyancy force in the formally unstable regions) be- 
comes zero in the surrounding stable regions, then adopts 
a suitable temperature stratification in the overshooting 
regions, and finally assumes straight mixing over-there. 
Since the Bressan et al. (1981) formalism makes use of 
the MLT, it expresses the mean free path of the convec- 
tive elements as / = A c x H p where H p is the local pressure 
scale height. 

Although some criticism has been advanced by Renzinj 
(1987(J] -who erroneously concludes that the "miscella- 
neous" of local and non local formalisms leads to overes- 
timating the overshooting distance- we remind the reader 
that overshooting is simply a logical consequence of the in- 
ertia principle, so that neglecting its existence would not 
be physically sound. It is worth recalling that convective 



overshooting is quite common in nature ( Deardorff et al. 
19690and it has been demonstrated in a number of stud- 



ies, including numerical simulations ( Freytag et al., 1996|) , 
that the penetration depth of convective elements into a 
formally stable region represents a non-negligible fraction 
of the size of the unstable zone. 

In addition to this, there are a number of astrophysical 
situations in which the hypothesis of substantial convec- 
tive overshooting has been found to offer better and more 



elegant solutions than other explanations (sec Bcrtclli 
et al., t_986 ; Chiosi ct al., 1992). Among others, it suffices 



1 It is worth recalling here that models of intermediate mass 
stars calculated with the classical Schwarzschild criterion are 
known to develop during the core He-burning phase the so- 
called He-semi-convective instability, i.e. a region surrounding 
the fully convective core in which the condition of neutrality 
Vij = V.4 is maintained by suitably modifying the profile of 
chemical composition. This is the physical analog of what hap- 
pens in massive stars during the core H-burning phase, which 
develop the so-called H-semi-convective instability. More de- 
tails on the phys ical origins of both H- and He-semi-convection, 
can be found in Chiosi et al. (1992). Let it suffice here to men- 



tion that in intermediate-mass stars the effects of semiconvec- 
tion are negligible, if compared to the effects of convective over- 
shooting. Hereinafter this type of stellar models are referred to 
as the classical, semi-convective models. 



stars (Bcrtclli ct al., 1995; Chiosi et al., 1992) 



Despite this, it has been often argued that unresolved 
binaries could mimic the effects of convective overshooting 
and easily account for the low ratio of red giant to main 
sequence stars observed in the young LMC clusters. This 
hypothesis has been investigated by many authors, both in 
open clusters of the Milky Way (see Carraro et al., 1994]) 
and young clusters o f the LMC ( Chiosi ct al., 1989a , bt 
Vallenari et al., 1992|) with somewhat contrasting results. 

In their study, Testa et al. (1999) emphasize the role of 
unresolved binaries in solving the problem, reaching the 
conclusion that overshooting is not needed. They obtain 
agreement with the observational data only by introducing 
a fraction of about 30 % of binary stars in a population of 
~ 100 Myr, and using the classical semi-convective models 
of Domingucz ct al. (1999| ). In contrast, models computed 
with "enlarged" convective cores (i.e. cores in which con- 
vection extends beyond the Schwarzschild criterion to sim- 
ulate overshooting) , lead to worse fits of the observational 
data, especially when binary stars arc included. 



Unfortunately, Testa et al. (1999) adopt a wrong pro- 



cedure to correct the observational star counts for photo- 
metric completeness (see the section below), thus obtain- 
ing differential (A^ms, shortly indicated with DLF) and 
integrated luminosity functions (EA^mSj shortly referred 
to as ILF) for the main sequen ce stars of NGC 18 66 that 
are inconsistent. Furthermore, Testa et al. (1999 ) fail to 
make use of the ILF normalized to the number of evolved 
stars (SA^ms /A^pmSi hereinafter indicated as N-ILF) in- 
troduced by phiosi et al. (198~9a ), which has been proved 
to be the only way to effectively discriminate between the 
two different evolutionary schemes. It is worth reminding 
the reader that the N-ILF is by definition proportional to 
the lifetime ratio tn/tn e . 

In this paper, first we correct the analysis of the ob- 
servational data and second we show that the conclusions 
reached by Testa et al. (1999| ) entirely depend on the kind 
of diagnostic they have adopted. Exploring the relative 
effects of core overshoot and unresolved binary star fre- 
quency on intermediate age cluster CMD's will be a side 
product of our analysis. 

Starting from the same observational data used by 
Testa et al. (1999[ )p| and correctly applying the complete- 
ness correction, in this paper we shall discriminate -by 
means of simulated CMDs and LFs of NGC 1866- between 
classical semi-convective and overshooting models. After 
a brief description of the observational data (Sect. ||), in 
Sect. | we compare the stellar models calculated accord- 
ing to the two schemes, with particular attention to the 
values of critical masses and lifetime ratios. These latter 
are indeed the crucial point to test the validity of stel- 
lar models. In Sect. [| we describe the isochrones and the 
methods used to simulate CMDs and LFs. In Sect. || we 
compare the results of these simulations with their obser- 



2 All the photometric data used in this study were kindly 
provided by Testa et al. (1999) already reduced and calibrated. 
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Fig. 1. CMD for stars i n the adopted sam ple, after statis- 
tical de-contamination ( Testa et al., 1999| ) 



vational counterparts. In Section ^ we thoroughly discuss 
the reasons why our conclusions differ from those by Testa] 
et al. fll999p . Finally, in Sect. @ we summarize the main 
results of this study. 

2. Observational data 



The observational data ( Testa et al., 1999 ) have been ob- 
tained with the 2.2 m ESO telescope at La Silla, Chile, on 
1993 January 27. Two fields have been observed in B and 
V pass-bands: the first field is centered on the cluster, and 
the second one is located at ~ 6' from it. The second field 
is utilized for statistical field subtraction: the corrected 
and calibrated CMD of Testa et al. (lgg'g] ) is presented 
in Fig. [I], in which the main sequence stars are separated 
from giants by dashed lines. 

2.1. Correction for completeness 



In thei r analysis, [Testa ct al. (1999 ) divide the sample in 

six concentric rings centered on the cluster, and calculate 
the completeness factors, A* , ring by ring. Since A* is equal 
to the number ratio of artificial stars added and recovered 
in a given magnitude bin, they should be used to correct 
the DLF - i.e. the number of stars per magnitude bin - 
separa t ely for each ring i : 



This clearly produces the correct DLF (iV l ) corr , but over- 
estimates the ILF and N-ILF at fainter magnitude bins. 
This can be appreciated by looking at the first 6 columns 
of Table [j], where we present the original star counts and 
the luminosity functions used by Testa et al. (199E ) 3 . In 
particular, it can be noticed that, starting from the un- 
corrected DLF and ILF of columns 2 and 3, the same 
"average" A c factors (column 4) have been used to derive 
the corrected DLF shown in column 5 (i.e. (A^Ms)corr = 
-^ms/Ac), and the corrected ILF shown in column 6 (i.e. 
£(AWcorr = SiV M s/A c ). The corrected DLF and ILF 
obtained in such a way are no longer consistent with each 
other. 

Noticing this, we decided to completely re-derive the 
luminosity functions for NGC 1866, starting from the sta- 



tistically decontaminated photometric data of Testa et al 
(1999| ) (see Fig. [I]), and the completeness factors presented 



in their Fig. 6. In short, we have obtained the DLF for each 
annulus around NGC 1866, and corrected them using the 
corresponding completeness factors separately, according 
to Eq. ([!]). Then, the DLFs for annulus 2 to 5 were added 



(N MS ) 



MS J COIT 



2=2 



MS7corr 



(2) 



producing the total DLF presented in column 8 of Table 
0. (For the sake of comparison, column 7 of the same table 
presents the counts obtained by adding stars from annu- 
lus 2 to 5, without applying the completeness correction.) 
^From this corrected DLF we derive the total ILF of col- 
umn 9 of the same table, by simply summing up the num- 
ber of stars above each magnitude bin. As can be readily 
seen, our final numbers for the DLF are roughly consis- 
tent with those of Testa ct al. (1999] ), but the two ILFs 
are sizeably different, especially for the faintest magnitude 
bins. 

There is a final remark to be made for the sake of clar- 
ity. Comparing columns (2) and (7), which contain the 
rough counts in the various magnitude bins, a marginal 
disagreement is evident: our counts do not exactly coin- 
cide with those by Testa et al. (1999). The difference is 



due to the slightly different criteria adopted to define the 



numbers Nms of stars per magnitude interval. Testa et al 
(1999|) follow a complicated scheme in which (i) the mean 



location of the MS together with its color dispersion a 
are derived as a function of the magnitude, and (ii) in ev- 
ery magnitude bin only stars falling within 7a from the 
mean MS are taken into account. In contrast, we prefer 
to consider all stars falling within the considered magni- 



tudc interval. As a consequence of this, the Nms of Testa 
ct al. (1999] ) are slightly smaller than ours because stars 



Nyxi 



where N z and (A^) corr stand for the observed and cor- 
rected inuinbers, respectively. 



(1) far away from the mean MS are neglected. 

In Fig. [2] we compare our results for the corrected 



Ho 



■wT 



ever, a s 

1999|) : they use the same A c factors to correct both the 
number of stars and the ILF at each magnitude bin 



ILF (continuous line) with those obtained by Testa et al 
(1999|) (dashed line). It is important to stress here that 



obstante error is made by ^sU el al.| Qnc ^ f rcasons why Qur results differ from thoge by 



PHrhese data have been kindly provided by V. Testa. 
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Table 1. Luminosity functions for main sequence stars in NGC 1866 as adopted in Testa et al. (1999): differential and integrated 
luminosity functions before the completeness correction (columns 2 ad 3), the completeness factors (column 4), differential and 
integrated luminosity functions after the completeness correction (column 5 and 6 ). The last three co lumns present the corrected 
luminosity functions as obtained independently by us, starting from the original Testa et al. (1999) data. 
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stars 
plied. 

to iV P 



et al. (1999), resides in the different (and erroneous) way 
of app] ying the completeness corrections used by them. 
Finally, we summarize in Table 2 the total number of 
ssquence (ETVms) and post main sequence (iVpMs) 
before and after the completeness correction is ap- 
The number of post main sequence stars amounts 
= 100. 



Table 2 

main 
correction 



seque nee 



Star counts for main sequence (ETVms) and post 
stars (,/Vpms), before and after completeness 



MS 



EiVMS A^PMS 



observed 
corrected 



4395 
4779 



100 
100 



3. Stellar tracks 

We have calculated two sets of stellar tracks with the clas- 
sical semi-convective scheme and initial chemical composi- 
tion [F=0.240,Z=0.004] and [y=0.250,Z=0.008], and one 
set of models with convective overshooting and composi- 
tion [Y=0. 240, Z=0.004]. The initial masses of the models 
range from 1.0 to 8 M Q . The input physics is the same 
as in Girardi et al. (2000 ) and subsequently updated by 
Salasnich et al. (2000|) ," to whom the reader should refer 



(2000) calculated with overshooting and chemical compo- 
sition [y =0.250,-2=0.008]. 

3.1. Prescriptions for semi-convection and 
overshooting 

The algorithm dealing with semi-convection in classical 



for all details. For purposes of compariso n, whenever re- 
quired we also utilize the stellar models by Salasnich et al. 



model s strictly follows the recipes given by [Alongi et al 



(1993 ) to whom the reader should refer for all details. 
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Fig. 2. Integrated Luminosity Function of the main se- 
quence stars of NGC 1866 before (dotted lines) and after 
(thick line) correcting for photometric completeness. The 
dashed line shows the same luminosity function calculated 
by [Testa et al. (1999| ) 



Convective overshooting is based on the formulation 
by Bressan et al. (1981) and the more recent revision by 
Bcrtelli et al. (1985[ ) and [Alongi et al. (1993| ). Suffice it 
to recall here that the overshooting parameter A is chosen 
to be A c =0.5 for core convection and A o =0.7 for envelope 
convection 



Alongi et al, 1993) 



For the sake of clarity and to avoid misunderstand- 
ing, we remind the reader that A c =0.5 in the formalism 



of Bressan ct al. (1981) is equivalent to A c =0.25 in the 
description of Maeder & Meynct (1991). More precisely, 
while Bressan ct al. (1981) define the overshooting dis- 
tance as the path traveled by a convective element across 
the Schwarzschild b order (half distance beneat h and half 
above the border), Maeder fc Meynet (1991 ) define it 
only as the distance above the border. Furthermore, the 
choice of A c =0.5 we have adopted is the same as in previ- 
ous studies aimed at generating stellar models best suited 
to interpret observat ional data of star c lu sters. To men- 
tion a f ew we recall Alongi et al. (1993]), Bressan et al 
(1993|)] [Maeder fe Meynet (19911 ), [Meynet et al. (1994| ) 



Chiosi et al. (1992)7 |Chiosi (199E ) and references therein. 
Exploring the effects of core overshooting for other choices 
of A c is beyond the scope of this study. 



3.2. Critical masses Mu e F, M, 



up 



As long ago pointe d out by |Chiosi et al. (1989a| ) , see also 
Chiosi et al. (1992| ), in the presence of core overshooting 



e 

i — — , 



o 



[ v =0.250, Z = 0.008] 

class 

overs] 




2.5M, 



4.4 



4.2 



4 

log T 



3.8 



3.6 



err 



Fig. 3. The stellar tracks for the 2.5 e 4 M© stars with com- 
position [Y=0.250, Z=0.008]. The solid lines are the stellar 
models with convective overshooting, whereas the dotted lines 
are the models with the classical extension of the convective 



Table 3. Critical masses Mn e F e M up - in units of M G 
for classical semi-convective and overshooting models. 



AW 



classical 
overshooting 



2.2 
1.8 



2.5 
2.0 



Table 4. Lifetimes for classical models with chemical com- 
position [y=0.250,Z=0.008]. The lifetimes are in Gyr. 



M/Mq 


tn 


t-Uc 




2.5 


0.4094 


0.2436 


0.595 


3.5 


0.1769 


0.1293 


0.389 


4.0 


0.1271 


0.0688 


0.340 


5.0 


0.0775 


0.0219 


0.283 


6.0 


0.0524 


0.0126 


0.241 


7.0 


0.0385 


0.0079 


0.208 


8.0 


0.0299 


0.0056 


0.189 



and core C-ignition in highly degenerate material occur, 
MhcF and M up respectively, get smaller. Table |^ sum- 
marizes the results obtained for classical semi-convective 
and overshooting models with metallicity Z=0.004 and 
Z=0.008. The same result holds good for both values of 
Z. 



the minimum initial masses, below which core He-flash 
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Table 5. The same as in Table ^ but for the overshooting 
models. 



M/Mq 




tHe 




2.5 


0.5438 


0.1376 


0.253 


3.5 


0.2308 


0.0380 


0.165 


4.0 


0.1674 


0.0237 


0.141 


5.0 


0.0999 


0.0116 


0.116 


6.0 


0.0672 


0.0065 


0.096 


7.0 


0.0487 


0.0043 


0.090 


8.0 


0.0376 


0.0032 


0.085 



3.3. Lifetimes and lifetime ratios 

Because of the increased size of the convective cores, the 
central H-burning lifetime in in the overshooting models 
is 30 % longer than in c lassical ones, this effe ct depending 
on the parameter A c ( |Brcssan ct al., 1981 ). The "over- 



luminosity" caused by overshooting during the core El- 
burning phase (see Fig. ||) still remains during the He- 
burning phase (for stars with initial mass greater than 
Mn e -p): consequently, the lifetime of the He-burning phase 
(£ho) gets shorter by about 40 — 45 %. This, combined 
with the longer in, results in a strong decrease (~ 55 — 
60%) of the ratio t Hc /t H . Tables | and | summarize the 
lifetimes for some of the Z=0.008 models calculated with 
the classical and overshooting scheme, respectively. 

4. Isochrones and synthetic CMDs 

4.1. Isochrones 

For the above four groups of evolutionary tracks, we 
have derived the isochrones, adopting the algorithm of 



"equivalent evolutionary points" developed by Bcrtclli 
et al. (1994). The age ranges from about 0.05 to 10 Gyr. 
Isochrones are calculated at Alogi = 0.01 intervals: this 
means that any two consecutive isochrones differ by ~ 2 % 
in their ages. Along each isochrone we give: the initial and 
actual stellar masses, the logarithm of surface luminosity, 
the effective temperature and surface gravity, the absolute 
bolometric magnitude, and the absolute magnitudes in the 
UBVRIJHK pass-bands. Finally we list the indefinite in- 
tegral over the initial mass M o f the initial mass fun ction 
(IMF). For any other detail see |Girardi ct al. (2000| ). 

Theoretical luminosities and effective temperatures are 
translated into magnitudes and colours in the Johnson- 



Cousins system, using the conversions of Bcrtclli et al 



(1994| ), which provides colours and bolometric corrections 
as a function of the e ffective temperature and gravity. The 
Bcrtclli ct al. (1994 ) conversions are based on the library 



of synthetic spectra obtained from theoretical models of 
stellar atmospheres by Kurucz (1992 ). 

We only note here, as illustrated in Fig. |[ that an 
isochrone of the overshooting scheme runs at higher lumi- 
nosities as compared to a "classical" isochrone of the same 
age, that is, in the presence of overshooting isochrones of 



-2 





Z=0.008, L=2xl0° yr 

classical 

overshooting 



1.5 



Fig. 4. Isochrones of same age (200 Myr) and chemical 
composition, [Y=0.250,Z=0.008] calculated with classical 
(dashed line) and overshooting models (solid line) 



older ages are required to fit the observational data of a 
cluster. 



4.2. Synthetic CMDs 

The synthetic CMDs are constructed by means of a Monte 
Carlo algorithm, which randomly distributes stars along 
a given isochrone according to evolutionary lifetimes and 
IMF. The code takes into account the photometric errors 



affecti ng the observational data (derived from Testa et al 
1999|) tables), by adding an artificial dispersion to the 



theoretical magnitudes and colours, and generates numer- 
ical simulations of CMDs allowing for effects due to: age 
spread of the cluster population, different slopes of the 
IMF, and presence of a certain fraction of unresolved bi- 
nary stars with mass ratio ranging in a given interval. 

Since any plausible IMF preferentially populates the 
main sequence region of an isochrone, in order to repro- 
duce the right proportion of main sequence to post main 
sequence stars -according to the lifetimes in the various 
evolutionary phases- we populate our synthetic CMDs 
until an assigned number of post main sequence stars 
^Vpms is matched: in particular in each simulation pre- 
sented here, we impose that the total number of red giant 
stars present in the data is matched. In such a way Npms is 
taken as the "normalization" parameter of the simulations 



(see Chiosi et al. 
|et al', 1991| ). 



1989a; Vallenari et al., 1991; Lattanzic 
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4.3. Luminosity functions 

Special care is paid to the integrated number of stars from 
the tip of main sequence band down to the current magni- 
tude interval, normalized to the number of red giant stars 
(i.e. the N-ILF). As pointed out by |Chiosi ct al. (1989a| ), 
the N-ILF can be directly compared with the ratio of life- 
times in the core H- and He-burning phases, in and tji e , 
respectively. At fixed -/Vpms, the ratio in each simulation 
EJVms/Wpms depends on the age -an older population 
clearly has a larger A^pms - and the theoretical model in 
use. As already noted in previous sections, the size of the 
convective cores strongly affects the ratio tn c /tn, and in 
turn, the ratio SJVms/^Vpms- Thus the N-ILF, when com- 
pared to the observational counterpart, reliably discrimi- 
nates classical from overshooting models. 

Concluding this section, it is worth remarking that re- 
producing synthetic CMDs with an assigned number of 
post main sequence stars is by far preferable to simulating 
CMDs with a fixed total number of stars brighter than a 
certain magnitude. As a matter of fact, the distribution of 
stars along the main sequence is driven almost exclusively 
by the IMF, evolutionary effects playing only a marginal 
role. Therefore, the LF (either DLF or ILF) of the main 
sequence stars, if not normalized to the number of giants, 
would only reflect the underlying IMF, and would not be 
affected by the effects of overshooting we want to test. 



5. The simulations 

In this section we present the CMDs and LFs resulting 
from our simulations, and compare them with the obser- 
vational ones, both for classical semi-convective and over- 
shooting models. 



[Y"=0.280,Z=0.020] in the present simulations we assume 
the composition [y=0.250,Z=0.008]. 



5.1.3. Distance modulus 



The distance modulus to the LMC is a hotly debated is- 
sue: we summarize here some of the most recent determi- 
nations of this parameter. Cepheid stars with Hipparcos 
parallaxes indicate (to - M) ~ 18.6 ( [Feast, 2000| ). The 
analysis of the expanding ring around SN1987A yields 
(m-M)o = 18.58 ±0.05 (|Panagia, 19981 ) ■ Th e simultane- 
ous stucl^ofJhc_Ccpheids stars and CMD of NGC 1866 
by [Bertelli et al. (1993Q yields (m-M) = 18.51 ±0.21 In 
a recent study on young LMC clusters ? utilize the value 
(to - M) = 18.45. We adopt here (m - M) = 18.5, 
as the "classical" value for the LMC distance modulus 
( Westerlund, 1997). Passing from (to 



M) = 18.5 to 
(to — M)q — 18.6 would not significantly change the re- 
sults of the present analysis. 



5.1.4. Reddening 

The study of the dust distribution across the LMC by 
means of photometric and spectroscopic data gives red- 
dening maps from which one derives an estimate of the 
mean value Eb-v = 0.16 ( pestrcichcr fc Schmidt-Kalcr 



1996), with a maximum value of 0.29 (reached in the re- 



gion of 30 Dor) and a minimum of 0.06 (Zaritsky et al 
1997|). While [Testa et al. (1999|) assume E B -v in the 



range between 0.06 to 0.10, we prefer to adopt the value 
0.10 which, t aking the total visual ext inction Ay to be 
3.1 x Eb-v ( Rieke fc Lebofsky, 1985 ), corresponds to 
Ay = 0.31. 



5.1. Adopted parameters 

Before presenting the simulations, we briefly summarize 
the values adopted for parameters such as reddening, 
metallicity, and distance modulus. 



5.1.1. IMF 



All the simulations are calculated adopting the galpeter 



(1955| ) law dN cx M~ a dM for the IMF and assuming, 



unless otherwise specified, the parameter a — 2.35. 



5.1.2. Metallicity 

Rough estimates of the metallicity of NGC 1866 de- 
rived from the BVI C colors of the Cepheids ( paldwell fc 



Coulson, 1985; Feast, 1989) yield [Fe/H]= -0.1 ±0.3 (see 



also Bcrtclli ct al., 1993| l, finally the recent analysis by 
Testa ct al. (1999b yielded [Fe/H]~ -0.35. This roughly 



corresponds to Z = 0.009. Therefore, improving upon the 
met allicity used in the old studies by Chiosi et al. (1989a ) 
and Brocato ct al. (1989] ) who adopted isochrones with 



5.2. Classical semi-convective models 

In Fig. [| is shown a synthetic CMD for the age of 76 Myr, 
and turn-off mass Mto = 5Mq. The termination magni- 
tude of the main sequence band fits the observational data. 
However, the luminosity interval spanned by the blue loop 
does not agree with the observed CMD. Theoretical lumi- 
nosities are indeed brighter than observed. Contrary- wise, 
if we impose that the mean luminosity of red giants is 
matched, the age increases to 110 Myr but the obser- 
vational termination magnitude of the main sequence is 
about 1 mag brighter than the theoretical one. In Fig. 
U we present the comparison between the observational 
(dashed line) and the theoretical N-ILF (solid line) for 
the population of 76 Myr in Fig. ||. The stellar models 
in use yield a ratio SA^ms/^pms lower than the obser- 
vational one. Since A^pms is taken equal to the number 
of red giants detected in NGC 1866, this implies that the 
classical semi-convective models produce more giants than 
observed (a long known result). 
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Fig. 5 



Fig. 7 



Fig. 5. Synthetic CMD constructed with classical models. The 
age is 76 Myr, the chemical composition is [y=0.250,Z=0.008]. 
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Fig. 6. Observational (dashed line) N-ILF — with the pois- 
sonian error bars — and theoretical N-ILF (thick line) for the 
case of classical models shown in Fig. ^. The age is of 76 Myr 
and the chemical composition is [Y—Q. 250,^=0.008] 



Fig. 7. Synthetic CMD for the age of 160 Myr, obtained 
from models with overshooting and chemical composition 
[Y=0.250,Z=0.i 
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Fig. 8. Observational (dashed line) N-ILF — with the pois- 
sonian error bars — and theoretical N-ILF (thick line) for the 
case of 160 Myr shown in Fig. [?] 



5.3. Models with overshooting 

As already anticipated commenting Fig. |J, the longer life- 
time and brighter luminosities of stellar models with over- 
shooting impose older isochrones to be adopted in the 
CMD simulations. Figs. [?] and || illustrate the synthetic 
CMD and the N-ILF, respectively, and compare them with 



their observational counterparts for a population 160 Myr 
old and turn-off mass of Mto = 4.1M©. Although a larger 
ratio EA^ms/^pms is obtained as compared with the clas- 
sical case, still the agreement between theory and obser- 
vations is not fully satisfactory: indeed, while the region of 
the blue loop is well fitted, the main sequence termination 
magnitude is about 1 mag fainter than observed. 
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Alternatively, if a younger age (~ 100 Myr) is adopted 
to match the main sequence termination luminosity, the- 
ory and observations would disagree in the morphology of 
the blue loop. Furthermore, it would predict too a high 
number of main sequence stars (ETVms — 10000) com- 
pared to the observed number of red giants (iVpMS = 100), 
i.e. the ratio SA^ms/^pms would be too high. 



5.4. Including binary stars 

The effect of unresolved binary stars (UBS) on the CMD 
of stellar cl usters has been inv e stigated in a number of 
studies (see Chiosi et al., 1989a ; Vallenari et al., 1991 ). 

The main effect of UBS on the main sequence mor- 
phology would be the appearance of a second sequence, 
running parallel to the main sequence, but systematically 
cooler and with brighter luminosities (~ 0.7 mag for bi- 
nary stars of equal mass). It goes without saying that in 
order to detect this secondary sequence the number of 
UBS must be significant. 

The percentage of UBS in open clusters of the Milky 
Way seems to amount to about 30 % - 50 % of the clus - 
ter p opulation ( |Mermilliod fc Mayor, 1989| ; |Carraro et al 



1994 ). Likely the same percentage of UBS could exist also 
in young clusters of the LMC. Strong observational evi- 



dence of this comes from the recent study of Elson et al 



(1998|) on NGC 1818, a young cluster similar to NGC 1866. 
In the CMD obtained with HST data, the existence of a 



"double" main sequence is soon evident. Elson et al. (1998 ) 
explain this feature by supposing the existence of a popu- 
lation of UBS amounting to 30-35 % of the total and with 
mass ratios in the range 0.7 to 1. Systems with mass ratios 
different from these cannot be excluded. In any case they 
would be hardly distinguishable from single objects. 

Basing on the study by Elson ct al. (1998), we include 



in our simulations the same fraction of binary stars, to 
understand how they alter the above conclusions both for 
classical and overshooting models. Unless otherwise spec- 
ified, the mass ratio q of the UBS included in the simula- 
tions falls in the i nterval 0.7 < q < 1 in agreement with 
Elson et al. (1998|) F|. 



The presence of UBS, shifting the main sequence ter- 
mination toward brighter luminosities, requires older ages. 
At the same time, the older the isochrone, the fainter is the 
blue loop. Therefore, under the action of the two effects, 
the luminosity gap between the main sequence termina- 
tion and red giants is reduced and the overall morphology 
of the CMD is better reproduced. 

5.4.1. Binary stars in the classical picture 

In Fig. H we present the synthetic CMD for the age of 105 
Myr (turn-off mass M TO = 4.3 M Q ), in which 30 % of UBS 



4 Detailed simulations of the CMD of N GC 1818, using the 
prescription for the population of UBS by Elson et al. (1998 ) 



can be found in Barmina (2001 ). They are not shown here for 
the sake of brevity. 



Fig. 9 



Fig. 9. Synthetic CMD constructed with classical models, 

chemical composition [Y=0.250,Z=0.008], and presence of 

UBS. The age is 105 Myr; the percentage of UBS amounts 
to 30 % 
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Fig. 10. Observational N-ILF (dashed line) — with the poisso- 
nian error bars — and theoretical ILF (thick line) for the case 
shown in Fig. ^[ classical models, age 105 Myr, and 30 % of 
UBS. 



are included. Compared to the results obtained simulating 
only single stars, the agreement with the observed CMD 
has much improved, but the theoretical N-ILF (the thick 
line in Fig. |Io|), clearly suffers the same difficulty already 
encountered neglecting the effects of binary stars, i.e. too 
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Fig. 11. Comparison between the observational N-ILF 
(dashed line) and the theoretical ones (thick lines) obtained in 
the classical picture for a population of 105 Myr and three dif- 
ferent slopes for the Salpeter IMF, namely a — 3, a — 3.5, a = 
4. The initial chemical composition, [V=0. 250,^=0. 008], the 
number of post main sequence stars (./Vpms = 100), and the 
fraction of binaries present in the sample (30 %) are constant 
in the three simulations 



low a ratio T,Nms/Npms- The same simulation has been 
repeated adopting the same percentage of binaries (30 %) 
but supposing that the mass ratio is distributed according 
to a Gaussian curve, as suggested by Testa et al. (1999| ). 
No significant difference is noticed both in the CMD mor- 
phology and N-ILF. 



5.4.2. Initial mass function 

How do the above results depend on the particular slope 
(a = 2.35 ) adopted for the Salpeter IMF? To cast light on 
this topic we have repeated our simulations varying a and 
keeping constant all other parameters. The analysis shows 
that while the CMD morphology is scarcely affected, the 
N-ILF is sensitive to variations of a. As illustrated in Fig. 
pr| , the N-ILF gets steeper at increasing a. We note that, in 
order to reach agreement between observational and theo- 
retical N-ILFs using the classical semi-convective models, 
one should adopt ev=(3.5 - 4), This would favor the forma- 
tion of low mass stars still located on the main sequence. 
The problem is whether or not such high values of a are 
acceptable. The present results re-confirm what already 



pointed out in the earlier studies by Becker & Mathews 
(1983) ) and phiosi et al. (198&J ). 



Fig. 12 



Fig. 12. Synthetic 
of 148 Myr 30% 
[Y=0.250,Z=0.i 
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CMD for overshooting models: age 
of UBS, and chemical composition 
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Fig. 13. Observational N-ILF (dashed line) -with the 
Poissonian error bars- and theoretical N-ILF (thick line) for 
the 148 Myr old population shown in Fig. |l2[ overshooting 
models and 30 % of UBS 



5.4.3. Binary stars in the overshooting picture 

The synthetic CMD presented in Fig. [l2] refers to stel- 
lar models with convective overshooting and simulates a 
population 148 Myr old (turn-off mass M T o = 4.2 M ) 
containing 30 % of UBS. We note that the luminosity gap 
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between the red giant region and the main sequence ter- 
mination has now disappeared. The CMD of NGC 1866 is 
perfectly reproduced. The most important result is, how- 
ever, the N-ILF shown in Fig. [l| which fully agrees with 
the observations. In other words, the ratio S7Vms/-^pms 
predicted by stellar models with convective overshoot to- 
gether with the presence of a certain fraction of binary 
stars fully agrees with that obtained from star counts. 



5.5. Effects of stochastic fluctuations in the IMF 

The discussion of the ILF we have presented so far, refers 
to a particular simulation of the CMD, which does not 
explicitly take into account possible effects of stochastic 
nature caused by the finite number of stars in a cluster. It 
is worth recalling here that the concept itself of IMF allows 
for important fluctuations to occur around the ideal distri- 
bution of stars as a function of the mass. The problem may 
be particularly severe for -/Vpms which in most cases is a 
small number. In other words the question to be addressed 
is whether stochastic effects may change iVpMS with re- 
spect to £-/Vms or conversely, whether iVpMS which has 
been used as the key constraint to our simulated CMDs 
and ILFs is generated by a unique population of main 
sequence stars or if many others are equally eligible. 

To cast light on this issue, we have performed many 
simulations of our reference CMD keeping constant all pa- 
rameters and letting the stochastic nature of the IMF 
develop thanks to the Monte Carlo technique we have 
adopted. To this aim, 100 simulations of the same CMD 
are calculated in order to get the mean N-ILF and its 
standard deviation a. 

For the sake of brevity, we illustrate here the results 
only for the cases of classical models (Fig. fl4|a) and over- 
shooting models (Fig. |l4]b) , in which 30 % of binary stars 
are included. No variations to the previous conclusions 
seem to appear. With classical semi-convective models the 
observed N-ILF (dashed line) runs above the ±2cr inter- 
val (dotted lines) around the mean N-ILF (thick line). In 
contrast, with the overshooting models theoretical and ob- 
servational N-ILFs agree at the ±2cr level. Considering the 
errors affecting the observational N-ILF, the agreement is 
extremely good. 



5.6. Age dispersion 

There is a final effect to be considered, i.e. the possibility 
that stars in the cluster are born over a finite time interval, 
i.e. that a significant age spread exists. 

Starting from our best-fit CMD, obtained from over- 
shooting models with age of ~ 148 Myr and including 
30% of binary stars, adding an age spread of about 20 
Myr (about 10%) would further improve the CMD, espe- 
cially in the region of red giant stars. This final CMD is 
presented in Fig. [l5|b, and compared with that of NGC 
1866 (Fig. |l5|a). This small spread in the age while giving 



the final make-up of the CMD does not affect the N-ILF, 
and the picture outlined in Fig. |l4] b still holds. 

5.7. Major uncertainties 

Looking at CMD of Figs. ||b and ||a, we see that while 
overall agreement is reached especially as far as the mag- 
nitudes and colors of the main sequence stars, and the 
mean luminosity of the stars in the blue loop are con- 
cerned, there is still a marginal disagreement in the colour 
range spanned by the red giant stars, i.e. the extension of 
the blue loop. This, however, is known to depend on sev- 
eral factors: (i) the overshooting parameters A c assumed 
for the nuclear region (see Chiosi et al., 1989a) and the 



one adopted for the envelope, A e (see Alongi et al., 1991); 
(ii) the rate of the nuclear reaction 12 C(a,7) le O -which 
determines the abundances of 12 C and 16 at the end of 
central He-burning phase- the higher the rate, the bluer 



the loop (Bertelli et al., 1985). It is worth noting, however, 
that these uncertainties on the extension of the blue loop 
do not affect the star counts at the base of the ILF. 

On the observational side, apart from photometric er- 
rors, completeness correction, and fraction of UBS that 
have already been taken into account, one has to check 
whether the results could depend in a crucial way on the 
particular calibration adopted to convert luminosities and 
effective temperatures into magnitudes and colours of the 
UBVRI system. To clarify th is point, we have als o used 
the empirical conversions by Alonso et al. (199£), who 



present tabulations of colours and bolometric corrections 
as a function of the effective temperature and metallicity 
[Fc/H] of the stars. However, comparing isochrones calcu- 
lated with this transformation and those obtained with the 
[Bertelli et al. (19*94 ) no appreciable difference is noticed. 

Some uncertainty could be related to the particular 
choice for the percentage of UBS (30 %) and their mass ra- 
tios 0.7 < q < 1 we have adopted. How could a change in 
these parameters affect our conclusions? To cast light on 
the issue we have performed additional simulations with 
different percentages and mass ratios. Increasing the bi- 
nary fraction up to about 40%, and extending the mass 
ratio interval down to 0.5 — 0.6, the CMD and N-ILF are 
still in agreement with the cluster data. However, assum- 
ing an even larger percentage of binary stars (> 40—45 %), 
the CMD gets worse as the main sequence band spreads 
too much toward red colors. 



6. Why do we disagree with [Testa et al. (1999|) ? 



The results of our analysis wholly disagree with those by 
Testa et al. (1999), despite the fact that we are using 
the same observational data for NGC 1866. Their main 
conclusion is indeed that classical stellar models without 
overshooting fit the CMD and ILF of NGC 1866 at the 
age of 100 Myr provided that a population of binary stars 
amounting to about 30 % of the total is included. In con- 
trast our main conclusion is that even including the same 
percentage of binary stars only models with overshooting 
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simultaneously reproduce the CMD, the ILF, and the cor- 
rect ratio £iVMs/-^PMS- The age now is about 150 Myr. 

What is the cause of the opposite conclusions? There 
are several factors concurring to the final result: 

(i) The different procedure applied to correct star 



counts for photometric completeness. The ILF of Testa 



et al. ( L999) is over-estimated due to the the fact that in- 
stead of summing up the DLF corrected for completeness 
in every magnitude bin, they apply the correction proce- 
dure to every step of the ILF. 

(ii) The different normalization procedures of the sim- 
ulations. In the present study, the normalization param- 
eter of the simulated CMD and ILF is the total (ob- 
served) number of red giants iVpMS- It is worth recalling 
that the N-ILF is given by T,Nms/^pms- The advantage 
of this type of normali zation has already been explained 
by Chiosi et al. (1989a| ) and briefly summarized in previ- 
ous sections. In contrast Testa et al. (1999|) assume that 
each simulation is complete as soon as a given total num- 
ber of stars with magnitude brighter than My = 2.6 is 
reached. They justify this choice by saying that their main 
sequence is well defined and complete up to three mag- 
nitudes fainter than the turn-off, and assume the total 
number of stars observed in this magnitude interval as 
the normalization parameter of the simulations instead of 
the giant (post main sequence) stars. 

(iii) This choice of the normalization factor is, how- 
ever, not physically sound because a deeply observed, well 
populated main sequence does not tell much about the in- 
ner structure, evolutionary rates, lifetimes etc. of its stars. 
Indeed the distribution of the stars along the main se- 
quence almost exclusively depends on the IMF. Therefore, 
the choice of a "normalization parameter" insensitive to 
the size of the convective core simply precludes all chances 
of discriminating between the two evolutionary schemes 
under examination. Simulations of CMDs and ILFs based 
on this criterion may fail to reproduce the correct ratio 

SA^ms/^pms- 

(iv) In relation to this, is also the fact that the to- 
tal number of stars brighter than My = 2.6 on which 
the simulations are normalized is determined by the pro- 
cedure of correcting for photometric incompleteness. We 
have already argued that the Testa et al. (1999|) method is 
inconsistent, and overestimates the total number of stars 
in the main sequence. Since their simulated CMDs and 
ILF do actually depend on this parameter, there is no 
way of becoming aware of the internal inconsistency. 

Considering all the differences between the present and 
Testa et al. (1999) analyses, we should recognize that the 
contrasting results are mainly due to the different way 
of normalizing the ILF: we choose a normalization fac- 
tor (number of giants) that is sensitive to overshooting, 
whereas Testa et al.'s (1999) choice (stars brighter than 
My = 2.6) is essentially insensitive to this effect. In fact, 
no clear discrimination between classical and overshooting 
models is expected with the latter method. 

(v) To better understand this point of controversy we 
have repeated our simulations strictly following the pro- 



cedure described by Testa et al. (1999): each one contains 
the same total number of stars brighter than V = 22 as 
the observational data, i.e. 4879. The results are displayed 
in Fig. limited to the few cases we have already exam- 
ined. The layout of Fig. ^ is organized in vertical and 
horizontal groups of panels. Panels in the same vertical 
row correspond to the case as indicated by the top head- 
ing of the row. Panels in horizontal rows show different 
types of LFs as indicated: the DLF in the bottom, the 
total ILF (main sequence and giant stars) in the middle 
panels, and the ILF for the sole main sequence stars in 
the top panels. The results can be commented as follows: 

— The DLFs in the bottom horizontal row make evident 
the luminosity gap between the main sequence ter- 
mination and the giant stars, located at My ~ —2. 
On the observational side, the gap is almost invisible, 
whereas on the theoretical side, the inclusion of bina- 
ries wipes out the gap in both classical and overshoot- 
ing models. All models seem to reproduce the main 
sequence part of the DLF (My — 1) equally well. 

— The total ILFs, i.e. inclusive of the post main sequence 
stars, are shown in the central horizontal row. For the 
models without binaries, the luminosity interval be- 
tween the main sequence termination and the giant 
stars appears as a plateau in the ILF at My ~ —2. To 
the left of this plateau, the detailed shape of the ILF 
is determined mainly by the location of giant stars: 
a slight over luminosity of the giants (with respect 
to observations) causes an apparent excess of stars at 
My < — 2, which can be noticed in the cases of classical 
models with and without binaries, and overshooting 
models without binaries. In the case of overshooting 
models with binaries, giant stars are predicted at the 
right magnitude level, and no such discrepancy can be 
noticed in the ILF. 

But it is important to notice that the "excess" sug- 
gested by the ILF does not necessarily correspond to a 
real excess of evolved stars. This is evident if we look 
at the total number of giants in each simulation (see 
the labels at the top of the figure) : in the case of classi- 
cal models there is a real excess of giants (158 and 189 
depending on whether or not binaries are included), 
which has no observational counterpart. The situation 
is the opposite for models with convective overshoot 
and no binary stars. Now the simulated number of gi- 
ants is 48 (about half of what observed) . Finally, mod- 
els with overshooting and binaries predict the correct 
number of post main sequence stars (120, to be com- 
pared with the 100 observed). 

— The ILF of main sequence stars (top horizontal row) 
does not show substantial differences passing from one 
case to another. This simply reflects the constancy of 
the IMF. There are small differences at the turn-off 
level, however too small to be significant. 

^From the plots in Fig. there is a hint that mod- 
els with overshooting and binary stars are in slightly bet- 
ter position, because they fit well the observed LFs in all 
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Fig. 16. Summary of the four cases under consideration: the heading of each column lists the type of stellar model (classical or 
overshooting), the pe rcentage of unreso lved binaries, the age, and the number of giant stars. The three horizontal rows display 
the LFs as shown in Testa et al. (1999): the top row is the integrated LF for the sole main sequence stars in the simulation, 
the middle row is the integrated LF for all stars in the simulation, finally the bottom row is the DLF for the whole sample. In 
each panel the solid line is the theoretical LF or DLF, whereas the dashed line is the same but for the observational data. Each 
simulation contains a total number of stars brighter than V — 22 equal to 4879 



panels. However, the most indicative result is not derived 
from the shape of the different LFs, but comes from the 
expected number of giants in the simulations: models with 
overshooting and binary stars are the only ones to predict 
the observed number of giant stars. 



In other words, with the kind of plot used by Testa 
et al. {1999|) , one cannot pin down the correct value of 
the ratio E7Vms/-^pms- Finally, we like to call attention 
to the fact that even with the Testa et al. (19%) method, 



no str ong indication is found that classical stellar models 
are to be preferred. 



7. Conclusions 



The main goal of this study is to compare with obser- 
vational data the predictions for CMD and ILF obtained 
from two different types of stellar models: the classical 
ones in which the extension of convective zones is defined 
by the Schwarzschild criterion (whenever required semi- 
convection is also taken into account), and the models 
with oversho oting (calculated accordingly to the Bressan 



et al. (198l| ) formalism), which allow the penetration of 



convective elements into the surrounding formally stable 
layers. The ultimate goal is to cast light on which type of 
stellar models find better correspondence with the physi- 
cal structure of real stars. 
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The main result that bears very much on the compari- 
son with the observational data is the net decrease -about 
55 — 60 %- of the lifetime ratio £hcAh passing from clas- 
sical to overshooting models. 

The analysis begins with correcting in the proper way 
the star counts and the LFs to get the LF of main sequence 
stars normalized to the number of evolved stars (N-ILF). 
The N-ILF yields the ratio £-/Vms/-^pms of main sequence 
to giant stars, and directly measures the ratio £hoAh (see 



Chiosi et al., 1989a) 



Many simulations of the CMD and LF of NGC 1866 
have been performed, in which not only the effect of the 
different stellar models but also those given by the pres- 
ence of binary stars is taken into account. The results can 
be summarized as follows: 

— The observed ILF cannot be matched at all by the 
classical models, showing too large a ratio of evolved 
to main sequence stars. The reason for the discrep- 
ancy between our results and the conclusion reached 
by Testa et al. (1999 ) has already been discussed and 



will not be repeated here 

— The situation is the opposite for the case of models 
with overshooting. It is slightly worse neglecting unre- 
solved binaries and fully satisfactory with binaries. 

— The presence of unresolved binaries in the simulations 
both for classical and overshooting models, allows to 
simultaneously match the luminosity of the main se- 
quence termination and giant stars. 

— Binary stars act on the main sequence morphology in 
a way roughly mimicking the effect of an enlarged H- 
burning core, i.e. shifting the turn-off toward brighter 
luminosities and forcing us to use older isochrones to 
fit the cluster data. However, they do not alter the 
ratio tHe/^H- 

— If all this can be taken as observational evidence for 
more extended mixing in the cores of hydrogen burning 
stars, the stellar models with convective overshooting 
constitute the simplest solution. They are indeed able 
to simultaneously match the CMD and ILF of NGC 
1866, thus re-confirming what already found in pre- 
vious studies (Chiosi et al., 1989a| Jb[ |Lattanzio et al. 
1990; |Vallenari et al., 199l[ |l992^ 

— By adopting the intrinsic distance modulus (m — 
M)o=18.5 and the initial chemical composition 
[y=0.250,Z=0.008], our best-fit of the CMD and ILF 
of NGC 1866 is for an age in the range 140 to 160 Myr, 
turn-off mass Mro=4-2 M , 30 to 40 % of binary stars 
with mass ratio q lying in the interval 0.6 < q < 1. 

— Finally, it is expected that adopting slightly different 
values for the distance modulus will not change the 
above conclusions. Looking at the turn-off luminosity 
and its relation with the age, we estimate that A(m — 
M)o=±0.1 with respect to the adopted value (m — 
M) = 18.5 would change the age (and age range) by 
At ~ ±10 Myr. Small variations in the percentage 
of binaries and/or age dispersion would immediately 
follow. 



New observational data will certainly improve the ex- 
act values of age, metallicity, binary star percentage so 
that the present values ought to be considered as provi- 
sional estimates. Although the new generation telescopes 
could make it possible to reach the innermost regions of 
the cluster, scarcely or not at all affected by field star 
contamination, thus providing us with better CMDs and 
ILFs, it is unlikely that the present conclusions will change 
significantly. Real stars should indeed possess bigger con- 
vective cores than predicted by the classical theory of stel- 
lar evolution. 

The civil war strikes again ... 
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Fig. 14. Mean N-ILF (thick line) and related 2a levels of 
expectancy (dotted lines) resulting from 100 simulations 
of the same population including 30 % of UBS. The chem- 
ical composition is [Y=0. 250, Z=0.008]. With the classical 
scheme (top panel) and age of 105 Myr the observational 
N-ILF (dashed line) is steeper than the theoretical one; for 
the overshooting scheme (bottom panel) and age of 148 
Myr, the observational N-ILF lies inside the ±2tr range 
around the theoretical prediction 
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Fig. 15a 



Fig. 15b 



Fig. 15. Comparison between the CMD of NGC 1866 ( [Testtj 
et al., 3999) (top panel) and our best-fit case (bottom panel), 
characterized by stellar models with convective overshooting, 
age in the interval 140-160 Myr, 30 % of UBS. The initial chem- 
ical composition is [K=0.250,Z=0.008]. The intrinsic distance 
modulus is (m — M)o = 18.5. 
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